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Topological insulators are new states of quantum matter with metallic edge/surface states. In 
this paper, we pointed out that there exists a new type of particle-hole symmetry-protected topo¬ 
logical insulator - topological hierarchy insulator (THI), a composite topological state of a (parent) 
topological insulator and its defect-induced topological mid-gap states. A particular type of THI is 
topological fractal insulator, that is a THI with self-similar topological structure. In the end, we 
discuss the possible experimental realizations of THIs. 


Recently, a new research topic focuses on searching 
for topologically nontrivial phases protected by symme¬ 
tries. Topological band insulator (TBI) provides an ex¬ 
ample of symmetry-protected topological states, in which 
there exist metallic edge/surface states [I|, 0. Accord¬ 
ing to the characterization of ” ten-fold way” of different 
symmetries Q, there are two types of TBIs in two dimen¬ 
sions: Z-type TBI with integer quantum Hall effect 0 
and Z 2 -type TBI with quantum spin-Hall effect[l|, d-Q. 
Then, people found that with extending the topological 
classification of band structures to include considering 
crystal point group symmetries, there exists new type of 
TBI - topological crystalline insulators [1|. Those addi¬ 
tional symmetries lead to a non-trivial topology of bulk 
wave functions and gapless edge/surface states. 


For TBIs, owning to the ’’holographic feature”, their 
nontrivial topological properties can be detected by prob¬ 
ing topological defects. For example, in two dimen¬ 
sional (2D) Z-type TBIs, 7r-flux traps mid-gap zero en¬ 
ergy states (zero modes) [lol - [l^ . For this topological 
state with a superlattice of 7r-fluxes, due to the poly¬ 
gon rule (Each fermion gains an accumulated phase shift 
\4>\ = (n — 2)tt 12 encircling around a smallest n-polygon 
rule [13), the defect-induced mid-gap states can be re¬ 
garded as an emergent ” TBI” with nontrivial topological 
properties, including the nonzero Chern number and the 
gapless edge states [l4|. 


In this paper we point out that a new kind of topo¬ 
logical insulator - topological hierarchy insulator may ex¬ 
ist, in which the topological properties are protected by 
translational and (generalized) particle-hole symmetry. 
In this system, the non-topological lattice defect - va¬ 
cancy traps zero mode(l5j|. For the system with vacancy- 
superlattice, the ground state becomes topological hier¬ 
archy insulator (THI). To show the nontrivial topological 
properties of THI, we take the Haldane model on square 
lattice and that on honeycomb lattice as examples. We 
found that the vacancy-superlattice could induce topo¬ 
logical mid-gap states and the low energy physics of 
the localized modes become that of an emergent ”TBI” 
on the parent TBI. In particular, the so-called topologi- 
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cal fractal insulator appears as a special THI with self¬ 
similar topological structure. 

The Haldane model on square lattice: Our start¬ 
ing point is the (spinless) Haldane model on square 
latticefH, [13 , of which the Hamiltonian is 

-ffparent = -|- -Ci -|- h.C.^ 

i^A 

+ 1 - icYc, + h.c.^ 

ieB 

-^'12 + h.c.J 

iGA 

+ t'Y + h.C.'j ( 1 ) 

iGB 

where Ci is the annihilation operator of the fermions at 
the site i. A and B label the sublattices, t and t' are 
the nearest neighbor (NN) and the next nearest neighbor 
(NNN) hoppings, respectively. For the Hamiltonian in 
Eq.(l), there exists 7r-flux in each square plaquette and 
^-flux in each triangular lattice. In this paper we set the 
lattice spacing a to be unit. 

The Haldane model on square lattice is a TBI with 

quantum anomalous Hall (QAH) effect. There is a finite 

energy gap near points ki = (0, 0) and k 2 = (tt, 0) which 

is A/ = 8t'. The Chern number Cparent of (parent) hal- 

dane model is 1 and the quantized Hall conductivity is 

£ 

h 

In particular, we point out that the Hamiltonian in 
Eq.(l) has the particle-hole (PH) symmetry. Under PH 
transformation, we have H^parent = —T’^H^parent'P where 
V = TZ ■ K, is the PH transformation operator [l^ [l8l| . 
Here, TZ is an operator that leads to Ci O (—l)®ct and /C 
is the complex conjugate operator. As a result, each en¬ 
ergy level with positive energy E is paired with an energy 
level with negative energy —E. Because the Hamiltonian 
in Eq.(l) is on a bipartite lattice, the quantum levels of 
the system with a single vacancy become an odd integer 
number. As a result, there must exist an unpaired elec¬ 
tronic state when we remove a lattice-site to create a va¬ 
cancy. Due to PH symmetry, the corresponding unpaired 
electronic state must have exactly zero energy [ 13 . We 
denote the wave-function of the zero mode by ipQ {ri — R) 
where R is the position of the vacancy. The quantum 
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states of the fermionic zero mode around a vacancy can 
be formally described in terms of the fermion Fock states 
{|0), |1)}. Here, |0), |1) denote the empty state, and the 
occupied state, respectively. 

Topological hierarchy insulator: We next study the 
Haldane model with vacancy-superlattice (VSL). When 
there are two vacancies nearby, the inter-vacancy quan¬ 
tum tunneling effect occurs and the fermionic zero modes 
on two vacancies couple. The lattice constant of the VSL 
is denoted by L, which is the distance between two vacan¬ 
cies nearby. When the vacancies locate at the different 
sublattices, the energy splitting AE is finite {L/a is an 
odd number); when the vacancies locate at the same sub¬ 
lattice, the energy splitting vanishes, AE = 0 {L/a is an 
even number). For the case of ^ < L, we can consider 
each vacancy as an isolated ’’atom” with localized elec¬ 
tronic states and use the effective tight-binding model 
to describe these quantum states induced by the VSL. 
Now, we superpose the localized states to obtain the sets 
of Wannier wave functions ipoifi — R)- In general, the 
effective tight-binding model of the localized states in¬ 
duced by the VSL becomes 

Hi-vl = - ^ TRR,d^^dR, - ^ TRR,d}^dR, ( 2 ) 

{R,R') {{R,R')) 

~ X/ "^RR'^^R^R' + 

{{{R,R'))) 

where dR is the fermionic annihilation operator of a lo¬ 
calized state on a vacancy R. Trr/ {T'rr,, is the 

hopping parameter between NN (NNN, NNNN) sites R 
and R' and is determined by the energy splitting AErri , 
\Trr,\ = lAifl {\Ti,^,\ = \AE'\, |T"^,| = |AF;"|). In 
general, due to <C \Trr,\ , \Trri\, we only consider 

the NN hopping term and NNN hopping term. Because 
the parent Hamiltonian in Eq.(l) has PH symmetry, the 
effective tight-binding model of the localized states in¬ 
duced by the VSL Hi-vl in Eq.(I2|) also has PH symme- 
Hy, Hi-vl = —V^Hi-vlV- 

Because the energy splitting AErri only determines 
the particle’s hopping amplitude \Trri\, we need to set¬ 
tle down the phase of Trri. In particular, to preserve 
PH symmetry, the total phase around a plaquette with 
four vacancies can be either tt or 0. An approach is to 
count the total flux number inside the plaquette. For 
the mid-gap states induced by VSL, every four vacancies 
form an L X L square, of which the total flux number is 
just L^/2 and the total phase around a plaquette is ttL^. 
After considering the compact condition, we have O/tt- 
flux inside each plaquette if is an even/odd number. 
In addition, we also use a numerical approach to check 
above prediction of the flux number inside the plaquette. 

Based on the Haldane model, we study the properties 
of a PH symmetry-protected TBI with VSL and intro¬ 
duce the concept of topological hierarchy insulator. To 
make it clear, we take the Haldane model with square- 
VSL as an example, of which the lattice constant L is 3a. 
See the illustration in Fig.l.(a). From above discussion. 
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FIG. 1: (Color online) (a) The illustration of square-VSLs 
on square lattice; (b) The particle distribution of zero modes 
around vacancies; (c) The DOS of the Haldane model with 
square-VSL for L = 3a: the mid-gap states are induced by the 
square-VSL. The inset shows the DOS of the parent Haldane 
model; (d) The edge states of mid-gap states. The parameter 
is t'/t — 0.2. 


there exists tt-Aux in each square plaquette of VSL and 
^-flux in each triangular VSL. Fig.l.(b) is the configu¬ 
ration of zero modes around vacancies. From Fig.l.(b), 
one can see that quantum states induced by VSL can be 
regarded as a new generation of lattice model, of which 
the vacancies play the role of the ” atoms”. 

As a result, we obtain an effective Hamiltonian to de¬ 
scribe the vacancy-induced mid-gap states, that reads 

Hi-vl = —T ^h+Lx'^R + '^^h+Ly^R + (^) 

ReA 

+ (^R+lJr - idli^LydR + h.c.) 

RGB 

+ ^ id^R.+ Lx+LvdR + dR+Lx-LydR + ^-C.) 

iGA 

~ ^ id^B.+Lx+LvdR + d^R+Lx-LydR + ^-C-) 

iGB 

where T, T' are effective NN and NNN hopping parame¬ 
ters, respectively. A and B label the sublattices of VSL. 
By comparing Eq. m and Eq.Q, one can see the corre¬ 
sponding relationship, Ci < —5- dR, t i —T, t' i — > —T'. 
For the defect-induced quantum states, the energy spec¬ 
trum and DOS are obviously similar to those of the par¬ 
ent Hamiltonian in Eq.(l). 

In Fig.I.(c) we show the DOS of the Haldane model 
with square-VSL for t'/t = 0.2 and L = 3a. In Fig.l.(c), 
the mid-gap (MG) energy bands appear and the DOS of 
the MG states is really similar to that of the parent TBI 
(the inset in Fig.I.(c)). We obtain the hopping param¬ 
eters by fitting the DOS of the defect-induced quantum 
states as T = 0.22t and T' = 0.07t. 

In addition, we will show the MG states induced by 
VSL have similar topological properties to those of par- 
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ent TBI. By the numerical approach, we found that the 
Chern number of the MG states is Ci_vl = which 
opposites to that of parent TBI as Ci-vl = —^parent- 
To characterize the THI, we define the filling factor, 
V = Nf jN where N is the total number of lattice sites 
and Nf is the number of fermions. The composite system 
is an insulator at filling factor v = i ot: u = 5(1 — f^)- 
At half filling (or v = ^), the total Chern number of the 
system is zero due to Ctotai = Cparent + Gi-vl; at filling 
factor i(l — jj), the total Chern number of the system 

is 1 due to Gtotal — Gparent- 

To illustrate the topological properties of the Haldane 
model with square-VSL, we then study its edge states. 
On one hand, at the filling factor i(l — ■^), the sys¬ 
tem becomes a TBI with Chern number I. As a re¬ 
sult, when the system has periodic boundary condition 
along y-direction but open boundary condition along x- 
direction, there exist the gapless edge states. On the 
other hand, we consider the half filling case, of which 
the system has zero Chern number. However, the system 
still has nontrivial topological properties. When the par¬ 
ent topological insulator has periodic boundary condition 
along both x-direction and y-direction, while the VSL has 
periodic boundary condition along y-direction but open 
boundary condition along x-direction, there may exist 
gapless edge states. See the numerical results in Fig.l.(d) 
for the case of L = 3a. We can see that there indeed exist 
gapless edge states on the boundaries of VSL. 

In general, the Haldane model with a square-VSL of 
odd number has nontrivial topological properties and 
is different from traditional TBIs. As a result, we call 
this type of PH symmetry-protected TBI with topological 
MG states to be topological hierarchy insulator. In other 
words, the point of view THI is given by the following 
equation 

THI = PH symmetry-protected TBI 
-I- L X L-square-VSL 

= Parent TBI -|- topological MG states. (4) 

Topological fractal insulator: We then study a particu¬ 
lar type of THI - topological fractal insulator (TFI), that 
is a THI with infinite generations of self-similar L x L- 
square-VSLs. See the illustration in Fig.2.(a). 

To construct a TFI, firstly we consider a THI with 
self-similar MG states. Because we can tune T'/T (or 
AE' / AE) by changing the NNN hopping t' of the parent 
TBI or the VSL constant L, the MG states can be self¬ 
similar to the parents states for the case of T/T' = t/t'. 
After numerical calculations, we obtain the self-similar 
condition in Fig.2.(b). See the dots in Fig.2.(b). Accord¬ 
ing to the phase diagram of Fig.2.(b), smaller/bigger tjt' 
leads to bigger/smaller T/T'. Now, under the self-similar 
condition T/T' = t/t', the effective Hamiltonian that de¬ 
scribes the MG states becomes iLi_vL = ct iJparent where 
a = T/< (a < 1) is an energy-scaling ratio. 

Next, we regard the vacancy-induced MG states as a 
parent TBI and consider an additional VSL on it. So, 
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FIG. 2: (Color online) (a) An illustration of the VSLs for 
TFI; (b) The self-similar condition; (c) The DOS of TFI. The 
parameter is t'/t — 0.2. The insets show the DOSs of the 
mid-gap states. 

using the recursive relation, we can see that there exist 
additional MG states inside the energy gap of first gen¬ 
eration MG states (we call it generation-1 MG states). 
Consequently, we get a new generation of topological 
MG states (we call such MG states in the energy gap of 
MG states as generation-2 MG states). Using the same 
method, we can construct a THI with n-generation MG 
states and eventually a THI with infinite generations of 
self-similar MG states that is just a TFI. Then, under the 
self-similar condition T/T' = t/t', the effective Hamilto¬ 
nian of the generation-n MG states is given by 

pj r a" Lfparent, n is an even number, , , 

Tln-VL = i n Pt* ■ JJ u 1^1 

a n IS an odd number. 

From the DOS in Fig.2.(c), we can see that the DOS of 
the MG states is always self-similar. 

In addition, we discuss the Chern number of the TFI. 
The TFI is an insulator at filling factor = ^(1 — j^) 
where m is an integer number, m > 1. The Chern 
number of generation-m MG states is known to be 
Cm-VL = (—1)'". As a result, when the filling factor 
is V = ^{1 — j^), the total Chern number of the system 
is 

m —1 

Gtotal = Gparent + X! = [1 + (6) 

n—1 

For example, when the filling factor is ^(1 — Tf) (m — 3), 
the total Chern number of the system is 1 due to Gtotai = 

Gparent + Gl_VL + G 2 -VL. 

It is known that a self-similar object looks ’’roughly” 
the same on different scales and fractal is a particu¬ 
larly self-similar object that exhibits a repeating pattern 
displaying at every scale. Topological fractal insulator 
provides a unique example of topological matters with 
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self-similarity, in which the topological properties always 
look the same on different energy scales or different 

length scales 

Topological hierarchy insulators on honeycomb lattice: 
We finally study the THI based on the Haldane model 
on honeycomb lattice, of which the Hamiltonian is given 

by 

Hparent = ~t ^ ^ cjcj — t ^ ^ 6 c\cj — p, ^ ^ cjci (7) 
(hi) <(hi)) i 

where Ci represents the fermion annihilation operator at 
site i. t and t' are the nearest neighbor (NN) and the next 
nearest neighbor (NNN) hoppings, respectively, is 
a complex phase along the NNN link, and we set the di¬ 
rection of the positive phase is clockwise = f)- M 
is the chemical potential. The Haldane model on honey¬ 
comb lattice is also a Chern insulator with nonzero Chern 
number, Cparent = 1- 

Since the Haldane model of Eq. © also satisfies the 
PH symmetry [H], iJparent = -)P^.ffparent)P, a vacancy 
induces a fermionic zero mode. In the followings, we fo¬ 
cus on the Haldane model with honeycomb-VSL shown 
in Fig.3.(a), of which the distance between two vacan¬ 
cies with shortest distance L is 5a. The particle density 
distribution of fermionic zero modes induced by VSL is 
shown in Fig.3.(b). 



FIG. 3: (Color online) (a) The illustration of honeycomb-VSL 
on honeycomb lattice. The distance between two vacancies is 
5a; (b) The particle distribution of zero modes around va¬ 
cancies; (c) The DOS of the Haldane model with honeycomb- 
VSL. The insets show the DOSs of the mid-gap states; (d) The 
edge states of mid-gap states. The parameter is t'/t — 0.1. 

Using similar approach, we derive an effective Hamil¬ 
tonian to describe the MG states 

i7i_VL =-T ^ ^ 

(L./) {(LJ» I 

( 8 ) 

where T, T' are effective NN and NNN hopping param¬ 
eters, respectively. /, J denote the positions of a vacan¬ 
cies. By comparing Eq. © and Eq. {|8j) , one can see the 


corresponding relationship to be Ci i—> di, t i — > T, 
t' i -5- T', (j)ij — —(filj. 

The DOS of the THI on honeycomb lattice is shown in 
Fig.3.(c). There exist MG states induced by honeycomb- 
VSL, of which the DOS is also similar to that of the par¬ 
ent Hamiltonian in Eq.(I71). Next, we calculate the Ghern 
number of the THI on honeycomb lattice and find that 
the Chern number of the MG states Ci_vl is —I which 
also opposites to the parent THI, Cparent = — Ci_vl- The 
THI is an insulator at filling factor v = ^ or v = At 
half filling (or v = i), the total Chern number of the 
system is zero; at filling factor zz = ||, the total Chern 
number of the system is 1. In addition, we study the 
edge states of the MG states. See the numerical results 
in Fig.3.(d). For a THI with n-generation MG states, 
the total Chern number of the system at hlling factor 
^=i(l-f^)is[l + (-l)'‘-i]/2. 

Discussion and conclusion: In this paper, based on 
the Haldane model on square lattice and that on honey¬ 
comb lattice, we studied the Chern insulators with VSL 
and found the vacancy-induced MG states have nontriv¬ 
ial topological properties. Then, we discover new types 
of TBI - THIs and TFIs. All these topological insulators 
are protected by particle-hole symmetry. In particular, 
TFIs provide a unique example of topological matters 
with self-similarity, in which the topological properties al¬ 
ways look the same on different energy scales. As a result, 
these exotic properties of THIs and TFIs will deepen our 
understanding of symmetry-protected topological states. 
In the future, similar topological states (THIs or TFIs) 
can be explored on other lattices (for example, the trian¬ 
gle lattice or kagome lattice). 

Finally, we address the relevant experimental realiza¬ 
tions of the THI. In condensed matter physics, it is dif¬ 
ficult to realize the Haldane model. However, the Kane- 
Mele (KM) model can be regarded as two copies of the 
Haldane model such that the spin up electron and the 
spin down electron exhibit opposite Chern nmber[I|. Be¬ 
cause the KM model also has PH symmetry, THI can be 
designed based on it. Ref. proposed that the KM 
model can be realized in a silicene-based material with 
a honeycomb lattice. A possible vacancy in silicene is a 
missing Si-atom that can be regarded as a vacancy in the 
KM model. For the silicene with periodic pattern of miss¬ 
ing Si-atoms, the THI may be possible realized. On the 
other hand, experimental realizations of quantum many- 
body systems in optical lattices have led to a chance 
to simulate a variety of topological states. Recently, in 
Ref.fl^, the (spinless) Haldane model on honeycomb op¬ 
tical lattice has been simulated in the cold atoms. Be¬ 
cause the VSL on optical lattice may be achieved by holo¬ 
graphic method, THI may be realized on optical lattice 
by using timely cold-atom technology. 
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